DESIGN AND STRESS ANALYSIS OF LOW-NOISE ADJUSTED BEARING 
CONTACT SPIRAL BEVEL GEARS 

A. Fuentes (1) , F.L. Litvin^, B.R. Mullins (3) , R. Woods (3) , R.F. Handschuh (4) , 

and D.G. Lewicki (4 ^ 

'T) Polytechnic University of Cartagena, Dept, of Mechanical Engineering 
C/Doctor Fleming, s/n, 30202, Cartagena, SPAIN 

( 2 ) Gear Research Center, University of Illinois at Chicago 
842 West Taylor Street, 60607 Chicago, Illinois, USA 

Bell Helicopter Textron, Inc., 76101 Fort Worth, Texas, USA 

U) NASA Glenn Research Center, U.S. Army Research Laboratory 
21000 Brookpark Road, Mail Stop 23-3, 44135 Cleveland, Ohio, USA 

Abstract 

An integrated computerized approach for design and stress analysis of low-noise spiral bevel 
gear drives with adjusted bearing contact is proposed. The procedure of computations is an 
iterative process that requires four separate procedures and provide: (a) a parabolic function of 
transmission errors that is able to reduce the effect of errors of alignment on noise and vibration, 
and (b) reduction of the shift of bearing contact caused by misalignment. Application of finite 
element analysis enables us to determine the contact and bending stresses and investigate the 
formation of the bearing contact. The design of finite element models and boundary conditions 
is automated and does not require intermediate CAD computer programs for application of 
general purpose computer program for finite element analysis. 


1 Introduction 

Design and stress analysis of spiral bevel gears is an important topic of research that has 
been performed by many scientists [1 - 6]. Reduction of noise and stabilization of bearing 
contact of misaligned spiral bevel gear drives is still a very challenging topic of research 
although the manufacturing companies (The Gleason Works (USA), Klingelnberg-Oerlikon 
(Germany-Switzerland)) have developed analysis tools and outstanding equipment for 
manufacture of such gear drives. 

The problem of design of spiral bevel gears is that machine tool settings for manufac- 
turing of spiral bevel gears are not standardized and have to be specially determined for 
each set of parameters of design to guarantee the required quality of the gear drives. A 
new approach is proposed for the solution of this problem that is based on the following 
considerations: 

(1) The gear machine-tool settings are considered as given (adapted, for instance, from 
the Gleason summary). The to-be-determined pinion machine-tool settings have to 
provide the observation of assigned conditions of meshing and contact of the gear 
drive. 
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(2) Low noise of the gear drive is achieved by application of a predesigned parabolic 
function of transmission errors of a limited maximum value of 6-8 arcsec. A pre- 
designed parabolic function of transmission errors is able to absorb almost linear 
discontinuous functions of transmission errors caused by errors of alignment. Such 
transmission errors are the source of high noise and vibration. 

(3) The provided orientation of the bearing contact has to reduce its shift caused by 
the errors of alignment of the gear drive. 

(4) The design procedure developed is an iterative process based on simultaneous ap- 
plication of local synthesis and tooth contact analysis (TCA). The local synthesis 
provides assigned conditions of meshing and contact at the mean contact point of 
tangency of pinion and gear tooth surfaces. The TCA computer program is ap- 
plied for simulation of the conditions of meshing and contact for the entire meshing 
process. 

Finite element method (FEM) is used for stress analysis and investigation of formation 
of the bearing contact. A model of three pairs of teeth complemented with the boundary 
conditions is applied for the finite element analysis (FEA). The development of contact 
models is automated and does not require application of CAD computer programs. 

2 Basic Ideas of Developed Approach 

Local Synthesis. The mean contact point is chosen on the gear tooth surface (Fig. 
1). The length of the major axis of the instantaneous contact ellipse at the mean contact 
point, 2a, the direction of the tangent to the contact path on gear tooth surface, i] 2t 
and the derivative of the gear ratio function m\ 2 are chosen at point M. The gear ratio 
function is given by mi 2 = u^/cu^ where a/ 1 ) and u/ 2 ) are the angular velocities of the 
pinion and gear rotations. The local synthesis program determines the pinion machine- 
tool settings considering as known the gear machine-tool settings and parameters a, i] 2t 
and 77i j 2 [7], The program requires solution of ten equations for ten unknowns but six 
of the ten equations are represented in echelon form. The algorithm of local synthesis 
includes relations between principal curvatures and directions proposed in [7, 8, 9]. 

Tooth Contact Analysis (TCA). The algorithm of tooth contact analysis is based 
on conditions of continuous tangency of pinion and gear tooth surfaces [7]. The TCA 
computer program determines the function of transmission errors and the bearing contact 
obtained for each iteration when the input variable parameters a, rj 2l m' 12 of the respective 
iteration are applied. 




Figure 1 : Illustration of parameters 772 and a applied for local synthesis. 

The computational procedure is divided into four separately applied procedures that 
are performed as follows: 

Procedure 1 : This procedure is directed at obtaining the assigned orientation of the 

bearing contact and is performed as follows: 

(a) The local synthesis and TCA are applied simultaneously. The varied parameter is 
m! n and parameters a and t ? 2 are taken as constant. The orientation of 772 is initially 
chosen for a longitudinally oriented bearing contact. The errors of alignment are 
taken equal to zero. 

(b) Using the output of TCA it becomes possible to obtain numerically the path of 
contact on gear tooth surface E 2 and determine its projection L T on plane T that 
is tangent to E 2 at point M (Fig. 1 ). 

(c) The goal of the iterative process is to obtain L ^ as a straight line for the process of 
meshing that is investigated for one cycle of meshing given by -1 r/Ni < d> 1 < 

This goal is achieved by variation of m r 12 and the sought-for solution is obtained 
analytically as follows: 

(i) The numerically obtained projection is represented by a polynomial function 

y t {x t , m'$) = Po (™'i 2 ) (i) + Pi(m' l 2 ) {i) x t + ft(m' 12 ) (,) x] (1) 


(ii) Variation of m'^ (* = 1,2,3, ... ,n) in the iterative process based on simulta- 
neous application of local synthesis and TCA, enables us to obtain a path of 
contact when d 2 — 0 and becomes a straight line. Figure 2 shows various 
lines L^\ L^\ and the desired shape L^\ 



Figure 2: Projections of various paths of contact L r on tangent plane T. 

(iii) Solution of Eq. (1) for ^ 2 (^ 12 )^ = 0- where m! n is the varied parameter, is 
based on the secant method [10]. 

Procedure 2: As a result of Procedure 1, the path of contact L ^ is represented by a 

straight line. However, the output of the TCA for the obtained function of transmission 
errors, function Ar/^ ^($i) where —tt/Ni < < ir/Ni, is of unfavorable shape and 

magnitude. The goal of Procedure 2 is transformation of A(^ n) (0i) into a parabolic 
function and limitation of the magnitude of maximum transmission errors. This goal is 
achieved by application of modified roll for pinion generation. 

Modified roll means that during the pinion generation the angles of rotation of the 
pinion and the cradle of the generating machine, designated as 'ipi and r ip c \ , respectively, 
are related as follows 


4>i\ipci) = m u ib cl - 4 J Vci - 1 (2) 

where m lc is the first derivative of function at ip el = 0 that is obtained by 

application of local synthesis [7], The superscript in Eq. (2) indicates that the j-th 
iteration is considered. 

Transformation of the function of transmission errors A^^i) is achieved as follows: 


(i) The function of transmission errors A^^^i) ' s the output of TCA obtained at 
the n-th iteration of Procedure 1 and is represented numerically. We represent 
A<4 n) (0i) as a polynomial function of the third order designated as 


= a ^ ) + + o> 3 J '(Pi, 
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(3> 


The designation j — 1,2, means that an iterative process for modification 

of A</4 n) (d>i) is considered. Function A^ j) (^i) = is obtained at the final 

iteration of Procedure 1. 
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The goal of Procedure 2 is to transform the function of transmission errors and 
obtain 
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where A<I> is the chosen maximum level of transmission errors. 


(iii) The goals mentioned above are obtained by variation of coefficients and b ^ 1 : of 
the function of modified roll. The secant method is applied for this purpose wherein 
variations of and are performed separately. 

Procedure 3: The goal of Procedure 3 is to reduce the shift of the bearing contact 

caused by errors of alignment and this is achieved by the proper change of orientation of 
L t assigned initially in Procedure 1 Procedure 3 is performed as follows: 

(i) Computer programs developed for local synthesis and TCA are again applied simul- 
taneously, but the expected errors of alignment are simulated. 


(ii) The effect of all errors of alignment on the shift of the bearing contact is investigated 
separately. The sensitivity of the shift of L T is reduced by the proper choice of 
parameter % (Fig. 1), that is the varied parameter of local synthesis in Procedure 
3. 


(iii) Application of procedure 3 enables us to obtain the resultant orientation of Lt 
that differs from the longitudinally one. The optimal orientation of L T depends on 
the design parameters of the gear drive, particularly on the gear ratio and on the 
required range of errors of alignment. 

After completion of Procedures 1, 2, and 3, the obtained pinion machine-tool settings 
guarantee that the designed gear drive is indeed a low-noise gear drive with reduced 
sensitivity to errors of alignment. 

Procedure 4: The goal of Procedure 4 is the stress analysis and investigation of formation 

of the bearing contact (see details in Section 5). 


3 Derivation of Equations of Gear Tooth Surfaces 

We recall that the machine-tool settings for the gear are considered as given and the to- 
be-derived equations allow the gear tooth surfaces to be determined. The head-cutter is 
provided with blades that are rotated about the Z 3 axis of the head-cutter (Fig. 3) during 
the process of generation. Both sides of the gear tooth are generated simultaneously. 
The profiles of the blade consist of two parts (Fig. 3): (i) of a straight lines, and (ii) of 
a rounded edge formed by circular arcs connected to the previously mentioned straight 
lines. The blades by rotation about the Zg axis form the head-cutter generating surfaces. 



Figure 3. Blade and generating cones for gear generating tool: (a) illustration of head- 
cutter blade; (b) and (c) generating tool cones for concave and convex sides. 


Applied Coordinate Systems. Coordinate systems S m2 , S a2 , and S b2 are the fixed 
ones and they are rigidly connected to the cutting machine (Fig. 4). The movable coor- 



dinate systems S 2 and S C2 are rigidly connected to the gear and the cradle, respectively. 
Coordinate system S g is rigidly connected to the gear head-cutter. It is considered that 
the head-cutter is a cone, and the rotation of the head cutter about the Z 0 axis does not 
affect the process of generation. The head-cutter is mounted on the cradle and coordinate 
system S g is rigidly connected to the cradle coordinate system S C2 . The cradle and the 
gear perform related rotations about the S m2 axis and the Sb 2 axis, respectively. Angles 
ib C2 and tb 2 are related and represent the current angles of rotation of the cradle and the 
gear. The ratio of gear roll is designated as m 2ci and is determined as 


m 2c2 = 


dtfo dtpcj 
dt ' dt 



Zm r Zc i 



Figure 4: Coordinate systems applied for the pinion generation: (a) and (b) illustration of 
tool installment for generation of righ- and left-hand pinions; (c) illustration of corrections 
of machine-tool settings. 

The installment of the tool on the cradle is determined by parameters S r , and q 2 , that 


are called radial distance and basic cradle angle. Parameters AXg 2 , A E mi , AX Dl , and 
7m 2 represent the settings of the gear. Figures 4(a) and (b) show the installment of the 
head-cutter for right-hand and left-hand gears, respectively. 


Derivation of Gear Tooth Surfaces. The head-cutter generating surface is rep- 
resented in coordinate system S g by vector function r g (s g ,0 g ) where s g and 9 g are the 
surface parameters. 

The family of generating surfaces is represented in coordinate system S 2 rigidly con- 
nected to the gear by the matrix equation 

Y 2 {s g ,9 g ,xp 2 ) = M 2g (ip 2 ) r g(s g ,6 g ) ( 7 ) 

where i <l> 2 is the generalized parameter of motion. 

The equation of meshing is represented as 


and determined as [7, 9] 

f2g(Sg,0 g ,lp 2 ) = 0 

(8) 

(dr lx dr 2 \ dra 

V 00 g 0s g ) d'lp 2 

(9) 

or as [7, 8, 9] 
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Here N g (s g ,O g ) is the normal to the head-cutter surface represented in coordinate system 
S g and v g ° ‘ is the relative velocity represented in S g . 

Equations (7) and (8) determine the gear tooth surface by three related parameters. 


4 Derivation of Equations of Pinion Tooth Sur- 
faces 

The two sides of the pinion tooth surfaces are generated separately. The machine-tool 
settings applied for generation of each tooth side are determined separately by application 
of Procedures 1, 2, and 3 previously mentioned. Profile blades of pinion head-cutters are 
represented in Fig. 5. 

Applied Coordinate Systems. Coordinate systems applied for generation of pinion 
are similar to those represented in Fig. 4. Coordinate systems Si and S 2 are rotated about 
the Z mi axis and Z bl axis, respectively (see Fig. 4 for equivalent coordinate systems Z m2 
and Z b2 for gear generation), and their rotations are related by a polynomial function 
Vi ('0c, ) wherein modified roll is applied (see below). The ratio of instantaneous angular 
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Figure 5: Blades and generating cones for pinion generating tool with straight blades: (a) 
convex side blade; (b) convex side generating cone; (c) concave side blade; (d) concave side 
generating cone. 

velocities of the pinion and the cradle is defined as m lc (tp i('</vi)) = ^ ( 1 ) (^ci)/w (c) - The 
magnitude of m lc is called ratio of roll or velocity ratio. Parameters AX Dl , AX Dl , A E mi , 
and 7 mi are the basic machine tool settings for pinion generation. 

Derivation of Pinion Tooth Surfaces. The pinion head-cutter surface is repre- 
sented by vector function r p ( s p ,6 p ) (Fig. 5) where ( s p ,6 p ) are the surface parameters. 
The family of head-cutters is represented in coordinate system Si rigidly connected to the 
pinion by the matrix equation 

r i{sp, Opj'tpci) — ) r p( 5 p> @p) (IT) 

Unlike the case of gear generation, modified roll is applied for generation of the pinion, 
and function ipi(ip ci ) relates the angles of rotation of the pinion and the cradle of the 


pinion generating machine by a polynomial but not linear function [see Eq. (2)]. The 
equation of meshing is represented as 

fpi 9p , 'ipci ) = 0 (12) 

and is determined by application of approaches similar to those represented by Eqs. (9) 
or (10). Equations (11) and (12) determine the pinion tooth surfaces by three related 
parameters. 

5 Application of Finite Element Analysis 

Application of finite element analysis permits: (i) investigation of the bearing contact 
wherein multiple sets of teeth may be in contact under load simultaneously, and (ii) 
determination of contact and bending stresses. 

Application of finite element method [11] requires the development of the finite ele- 
ment model mesh, the definition of possible contacting surfaces, and the establishment 
of boundary conditions for loading the gear drive with the given torque. Finite element 
analysis is performed by application of general purpose computer program [12]. 

The development of the solid models and finite element meshes using CAD computer 
programs is expensive, requires skilled users of those computer programs, and has to be 
done for every case of gear geometry and position of meshing to be investigated. 

The authors have developed a modified approach to perform the finite element analysis 
that has the following advantages: 

(a) Computational procedures for synthesis, analysis, and generation of finite element 
models are integrated in developed computer programs. Graphical interpretation of 
the output is obtained by using commercially available graphical subroutines. 

(b) The generation of the finite element mesh required for FEA is performed automat- 
ically using the equations of the surfaces of the tooth and its portion of the rim. 
Nodes of the finite element mesh that belong to the tooth surfaces of pinion (gear) 
model guarantee to simulate points of the real tooth surfaces of the pinion (gear). 
Loss of accuracy due to the development of solid models using CAD computer 
programs is avoided. The boundary conditions for the pinion and the gear are set 
automatically as well. 

(c) The formation of the bearing contact for a cycle of meshing can be investigated 
and edge contact can be discovered and avoided. 

(d) Setting of boundary conditions for gear and pinion is obtained automatically. Nodes 
on the sides and bottom part of the rim portion of the gear are considered as fixed. 



Nodes on the two sides and bottom part of the rim portion of the pinion build a 
rigid surface. The rigid body reference node is located on the pinion axis of rotation 
with all degrees of freedom except the rotation around the axis of rotation of the 
pinion that is fixed to zero. The torque is applied directly to the remaining degree 
of freedom of the rigid body reference node. 

(e) Definition of contacting surfaces for the contact algorithm of the finite element 
computer program [12] is automatic as well and requires definition of the master 
and slave surfaces. Generally, the master surface is chosen as the surface of the 
stiffer body or as the surface with the coarser mesh if the two surfaces are located 
on structures with comparable stiffness. Finite element mesh for the pinion is more 
refined than the one for the gear due to larger curvatures. The gear and pinion 
tooth surfaces are considered as the master and slave surfaces, respectively, for the 
contact algorithm. 

6 Numerical Examples 

Three examples of design of face-milled spiral bevel gears with various design parameters 
have been investigated. The output of computations for the mentioned examples are 
represented in Fig. 6 and shows: 

(i) The bearing contact has a longitudinal direction in Example 1 for an aligned gear 
drive [Fig. 6(a)]. 

(ii) The bearing contact in Example 2 [Fig. 6(b)] is deviated from the longitudinal 
direction for reduction of the shift of bearing contact caused by errors of alignment. 

(iii) Example 3 of design [Fig. 6(c)] shows that the gear drive with a gear ratio close to 
1 {N 2 /Ni = 26/17) is very sensitive to errors of alignment. The shift of bearing 
contact due to misalignment is reduced since larger deviation of the bearing contact 
from the longitudinal direction has been provided. 

(iv) The function of transmission errors [Fig. 6(d)] for all three examples is indeed a 
parabolic function up to 8 arcsec of maximum amount. 

Finite element analysis has been performed for all three examples. Figure 7 shows the 
formation of the bearing contact at the heel position of the gear for one of the studied 
gear drives. 



(d) 


Figure 6: Bearing contact on gear tooth surfaces and predesigned function of transmission 
errors for: (a) longitudinally oriented bearing contact for gear drive of example 1; (b) 
adjusted bearing contact for gear drive of example 2; (c) adjusted bearing contact for’gear 
drive of example 3; (d) function of transmission errors for examples 1, 2, and 3. 

7 Conclusions 

Based on the study conducted the following conclusions can be made: 

(1) An integrated computerized approach has been developed for design of low-noise 
spiral bevel gears with an adjusted bearing contact based on the following ideas: 

(i) A parabolic function of transmission errors of limited value of maximum trans- 
mission errors is applied. Such a function is able to absorb linear discontinuous 
functions of transmission errors caused by misalignment. 
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Figure 7: Bearing contact at. heel contact point on gear tooth surfaces for spiral bevel gear 
drive. 

(ii) The orientation of the bearing contact is adjusted for reduction of the shift of 
the bearing contact caused by errors of alignment. 

(iii) The approach developed is an iterative procedure based on simultaneous ap- 
plication of local synthesis and tooth contact analysis (TCA) using modified 
roll for pinion generation. 

(iv) The contacting model for finite element analysis is formed by three teeth 
and the boundary conditions, it is automatically designed, and does not need 
intermediate CAD computer programs for application of finite element analysis. 

(v) The same computer language is applied for numerical computations performed 
for all stages of design, and automatic generation of finite element models. 

(2) The approach developed may be also applied for design of formate cut spiral bevel 

gears, hypoid gear drives and other types of gear drives. 
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